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ABSTRACT 

The hierarchical formation scenario for the stellar halo requires the accretion and disruption of dwarf galax- 
ies, yet low-metallicity halo stars are enriched in a-elements compared to similar, low-metallicity stars in 
dwarf spheroidal (dSph) galaxies. We address this primary challenge for the hierarchical formation scenario 
for the stellar halo by combining chemical evolution modelling with cosmologically-motivated mass accretion 
histories for the Milky Way dark halo and its satellites. We demonstrate that stellar halo and dwarf galaxy 
abundance patterns can be explained naturally within the ACDM framework. Our solution relies fundamen- 
tally on the ACDM model prediction that the majority of the stars in the stellar halo were formed within a 
few relatively massive, ~5x 1O 1O M , dwarf irregular (dlrr) - size dark matter halos, which were accreted and 
destroyed ~ 10 Gyr in the past. These systems necessarily have short-lived, rapid star formation histories, are 
enriched primarily by Type II supernovae, and host stars with enhanced [a/Fe] abundances. In contrast, dwarf 
spheroidal galaxies exist within low-mass dark matter hosts of ~ 10 9 M©, where supernovae winds are impor- 
tant in setting the intermediate [a/Fe] ratios observed. Our model includes enrichment from Type la and Type II 
supernovae as well as stellar winds, and includes a physically-motivated supernovae feedback prescription cal- 
ibrated to reproduce the local dwarf galaxy stellar mass - metallicity relation. We use representative examples 
of the type of dark matter halos we expect to host a destroyed "stellar halo progenitor" dwarf, a surviving dlrr, 
and a surviving dSph galaxy, and show that their derived abundance patterns, stellar masses, and gas masses 
are consistent with those observed for each type of system. Our model also self-consistently reproduces the 
observed stellar mass - v c i rc relation for local group satellites and produces the correct cumulative mass for the 
Milky Way stellar halo. We predict that the lowest metallicity stars in intermediate-mass dlrr galaxies like the 
SMC and LMC should follow abundance patterns similar to that observed in the stellar halo. Searches for ac- 
creted, disrupted, low-mass dwarfs may be enhanced by searching for unbound stars with dSph-like chemical 
abundance patterns. 

Subject headings: galaxies: formation - galaxies: evolution - Galaxy: halo - Galaxy: abundances - galaxies: 
dwarf 



1. INTRODUCTION 

While evidence in support of a ACDM universe contin- 
ues to mount fro m high-precision, large-scale cosmological 
observations (e . g. ISoergel et alj|2003|: |Blakeslee et ap l2003l 
iPen et alJl200l lLahav et al.ll2002HZehavi et alJl2002fc . a va- 
riety of data on galactic and subgal actic scales has rem ained 
more difficult to explain (e.g.lSimo n et alj|2005|: |Mo ore et al.l 
ll999HKlvpin et alJl999HNavarro & SteinmetJ 19971 and ref- 
erences therein). Among the biggest questions in cosmology 
is whether ACDM needs to be modified substantia lly in order 
to account for galactic-scale data (e.g. ISoergel & Steinhardtl 
l2000t lKaDlingh at et al l2000UColm et all2000l) or if physical 
and dynamical processes can expl ain the data in the context 
of the standard paradigm (e.g. Bullock et al. 2000; lAlam et akl 
120021 iRobertso n et al.ll2004l) . The resolution of this question 
requires testing our cosmological model against the full realm 
of existing subgalactic data. In this paper, we confront the rich 
data set of elemental abundances in galactic halo and Local 
Group dwarf galaxy stars with models set within the ACDM 
cosmology. We conclude that the chemical abundance pattern 
in dwarf irregulars, dwarf spheroidals, and the Milky Way 
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stellar halo can be explained naturally in the context of our 
hierarchical cosmology. 

One of the fundamental tenets of the CDM paradigm is that 
dark matter halos form hierarchically, via a series of mergers 
with smaller halos. This prediction gives rise to the natural 
expectation that the ste llar halo is formed from disrupted, ac 
creted dwarf galaxies (Qphnston et al. 1996; Helmi & Whit 
19991 iBullock et all 1200 It iBullock & Johnston! 120051 
Font et aljT27)05"ln^ reparatioril). Indeed, this pictu re is 



remarkably similar to the ideas of ISearle & Zinnl lll978l) who 
argued for "chaotic accretion" based purely on observational 
grounds. The discovery of fossil evidence of these accretion 
events via the identification of substructure within the stellar 
halo of our Galaxy may provide the only direct evidence 
that structure formation is hierarchical on small-scales 
iBullock & Johnston! 120051). Searche s for this substructure 
are u nder way (e.g. |Ne wberg_eLalJ [2QQ2t |M^rrisQ_n_eLaI] 
[200l [lKlniiian et al.l2004UBrown et alJ2004llMaiewski et alJ 
120041) . 

One of the main problems facing the hierarchical accretion 
picture for the stellar halo is that abundance patterns of sur- 
viving local group galax ies do not match those of s tars in 
the stellar halo (see e.g. IWvsdl2004 IVenn et al J 120041 
Stars in the stellar halo tend to be old, metal-poor, and a- 
enriched, with [a/Fe] ~ +0.3. The stellar halo abundance pat- 
tern contrasts with the younger stellar populations of local, 
massive dwarf irregular (dlrr) galaxies like the Large Mag- 
ellanic Cloud (LMC) and Small Magellanic Cloud (SMC), 
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whi ch are more m etal enriched and have nea r-solar fg/Fel ra- 
tios ( Jasniewicz & Thevenin 1994; Hill et al. 1995; Hill 1997; 
[ Hill et all 1997H Venn! 1999UCole et al.l2000UHill et alJ200ft 
iKorn et al.ll2002l ISmithl 120021) . Additionally, stars in dwarf 
spheroidal (dSph) galaxies around the Milky Way have inter- 
mediate a-element abundanc es at low-meta ll icitv, Ta/Fel ~ 
+0.15 (ISmecker-Hane & Mc William! 1 19991: I Shetrone etaT 
1 2001 1 : ISmecker-Hane & Mc William! [20031 Ishetrone et al 
l2003t iVenn et alJ 120041 iGeisler et alJ 12005b . The discrep- 
ancy between the ages and metal abundances of dwarfs 
and the stellar halo is often cited as evidence agai nst the 
hierarchical ACDM pr ediction (e.g . lUnavane et alJ 1199' 
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Venn et al. 2004), and has even been used to argue for a re- 
vision of the ACDM cosmology (Wvse 2004). The problem 
is clear: how can one build the stellar halo from destroyed 
dwarf galaxies when the stars in low-mass galaxies appear to 
be chemically distinct from the stars in the halo? 

Our solution to this conundrum relies fundamentally on two 
predictions that follow directly from the ACDM cosmology: 
1) most of the mass in the Milky Way halo was acquired via 
mergers with massive (~5x 1O 1O M ) "dlrr-type" dark matter 
halos; and 2) most of these events occurred quite early on, at a 
look-back time of ~ 10 Gyr. It follows naturally that the ma- 
jority of stars in the stellar hal o were formed in in termediate- 
mass halos of this type (Bullock & Johnston 2005). The asso- 
ciated galaxies were disrupted soon after accretion and before 
significant chemical evolution. 

Within the context of our scenario, the observed tendency 
for halo stars to be metal poor and a-enhanced compared to 
stars in surviving dlrr galaxies follows from standard assump- 
tions about Type la vs. Type II supernovae enrichment. Type 
II supernovae provide a-enhanced enrichment within ^10 
Myr of a star-formation episode (see |4}, while the iron-rich 
nucleosynthetic products from Type la supernovae contribute 
much later, at ~ lGyr. The massive, "stellar halo progeni- 
tor" halos are accreted early and destroyed soon after. These 
systems are forced to have short, truncated star formation his- 
tories, and naturally host stars that are predominantly Type-II 
enriched and a-enhanced. In contrast, surviving dlrr galax- 
ies are accreted much later than their destroyed counterparts 
(this is what allows them to survive, see ^2J, form stars for 
many billions of years, and have substantial enrichment from 
both Type la and Type II supernovae. Finally, lower-mass, 
dwarf spheroidal galaxies inhabit halos with much shallower 
potential wells (< 10 9 M Q ) than the massive dlrr galaxies. We 
argue below that blow-out feedback from supernova give rise 
to both the low-metallicities observed in these systems as well 
as their intermediate [a/Fe] ratios. 

In what follows, we develop a model to track the chem- 
ical enrichment histories in dwarf galaxies forming in cos- 
mologically self-consistent dark matter halos. We use three 
example mass accretion histories. Two are chosen to be rep- 
resentative of the kind of dark matter halos we expect to host 
surviving dSph and dlrr satellites, based on ACDM predic- 
tions. A third example accretion history is typical of the type 
that is accreted and destroyed early to contribute to the forma- 
tion of the halo. We track chemical enrichment by including 
contributions from Type la and Type II supernovae as well as 
stellar winds. Our star formation and feedback model self- 
consistently reproduces: 1) the total mass of the stellar halo 
(Bullock & Johnston 2005); 2) the stellar mass-v c i, c relation 
expected for Local G roup dwarfs; and 3) th e mass - metallic - 
ity relation found bv lDekel & Wool (120031) (see We go 



on to use this model to interpret the chemical abundance pat- 
terns of the stellar halo, dwarf spheroidals, and dwarf irregular 
galaxies dVenn et all2004l) . 

We note that there exists an extensiv e literature of 
galac tic chemical evolution mod elling ijyan de n Bergh 
19621 iLarsonl Il972t [Tinslevl Il972t iTalbot & Arnettl 119731: 
Pagel & Patchetll975tlHartwickll976tlTinslevll980tlTwarog| 
1980J IMatteucci & Greggiol I1986L iMatteucci & Tor nambe 
1987t IMatteucci & Francois! Il989t IWheeler et alJ J1989J: 
Ferrini et alll992tlTimrnes et alJl995t|Tsuiimoto et al.ll995t 
pevriendt et al.lll999t iFerreras & SiipboOOt IChiappini et all 
120031 lLanfranchi & Matteuccil 120041) . We aim neither to 
critique nor supersede this previous work. Our goal is simply 
to approach galactic chemical evolution with a hierarchically 
driven and cosmologically consistent method while retaining 
complexity in our chemical enrichment code (see also 
Brook et al. 2004a b). By placing the chemical enrichment of 
dwarf galaxies into the context of their expected cosmological 
formation and subsequent fate, we hope to gain some insight 
on interpreting the rich data sets being compiled from Local 
Group observations. 

In the next section, we provide an overview of the ex- 
pected formation histories of Milky Way-type dark matter ha- 
los, the accretion and disruption histories of their accreted 
subhalos, and the mass accretion histories of the satellite ha- 
los themselves. In ^3]we discuss our star formation law and 
present our three example mass accretion histories used to fa- 
cilitate our comparisons to dlrr, dSph, and stellar halo abun- 
dance data. We detail our stellar yield modelling in 5j4] and 
our chemical evolution and feedback calibration in 50 We 
present our results in ^6] and summarize and conclude in 50 
Throughout this paper we assume a flat ACDM cosmology 
with fi m = 0.3, Q, A = 0.7, h = 0.7, and cr 8 = 0.9. 

2. COSMOLOGICAL CONTEXT 

We first review several relevant expectations for the for- 
mation of galaxy-size dark matter h alos within the ACDM 
framework. The reader is referred to Wechsler et al. (2002); 
IZentner & Bullockl (2003) and lZentner et alJ J2005I) for more 
detailed discussions. 

2.1. Hierarchical growth of the Milky Way dark halo 

Consider a galaxy-size dark matter halo of mass Mh at z = 0. 
Mass growth in halos of this type proceeds by the merging and 
accretion of a wide mass spectrum of "subhalos" or "satel- 
lite halos". Integrated over the lifetime of the host halo, the 
accreted spectrum of satellite halo masses, M s , rises softly 

—1/2 

towards the low-mass end, dN/dM s oc M s 1 and falls off 
sharply above a characteristic cutoff mass ~ O.lMh. Owing 
to the shallow slope, most of the mass accreted onto the host 
is associated with subhalos of roughly the cutoff mass. 

The Milky Way's dark matter halo has a virial mass of 
Mh ~ 10 12 M Q at z = and is thus expected to have accreted 
most of its mass via the accr etion of in termedia te size sys- 
tems M s ~ (1-10) x 10 10 M Q (IZentner & Bullockl2003l) . Ac- 
creted subhalos lose mass as a result of violent interactions 
with the host halo system, and are typically destroyed within 
a few dynamical times. The stripped material incorporates 
into and effectively bu ilds the b ackground d ark ma tter halo of 
the host. The model of Bullock & Johnston ( 2005) that serves 
as our cosmological framework accounts for these dynamical 
effects, in addition to the redshift evolution of the characteris- 
tic mass and densities of merging systems. 
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The time line of accretion activity is sensitive to the host 
halo mass and the governing cosmological parameters. For 
Milky-Way size halos in the standard ACDM model, the ac- 
cretion rate of dwarf-galaxy halos peaks ~ 10 Gyr in the past, 
and slows down considerably at a look-ba ck time of ~ 8 Gyr, 
althou gh there is scatter from halo to halo (Zentner & Bullock 
120031) . Satellite halos that get accreted during the earliest 
epochs have time to be destro yed by i nteractions with the 
host galaxy's environment (see Zentner et al. 2005). The in- 
creased efficiency of dynamical friction with satellite mass 
causes more massive subhalos to lose mass more quickly than 
smaller subhalos. Most of the subhalos that survive until 
the present day were accreted during the slow-growth phase, 
within the past ~ 5 Gyr, and therefore represent a temporally 
biased population compared to their destroyed counterparts. 
We emphasize that this temporal bias between surviving and 
destroyed systems implies that the two populations have in- 
trinsically different formation histories: the destroyed popu- 
lation of halos are expected to have formed more rapidly than 
the subhalos that survive. 

We expect that the most massive, accreted satellite halos 
will host dwarf galaxies, and that most of these systems will 
be destroyed by tidal interactions. Just as the dissipationless 
dark halo of the host is built from the dark matter that was 
originally associated with accreted and destroyed subhalos, 
the stellar halo is built from the dissipationless stellar mat- 
ter that wa s associated with th e same population of destroyed 
galaxies ( Bullock et a Given this scenario for stel- 

lar halo formation, two important implications follow directly 
from the ACDM predictions discussed above: 

1. The majority of the stars contained in the Milky 
Way stellar halo originated in ~ 5 x 10 10 M Q "dwarf 
irregular-type" dark matter halos. 

2. The surviving satellite galaxies around the Milky Way 
represent a biased population: they typically were ac- 
creted later than the dwarf galaxies that were destroyed 
to make the halo. 

The second point implies that the mass accretion and star for- 
mation histories of surviving satellites are more gradual than 
their destroyed counterparts. 

2.2. Dark matter halo mass accretion histories 

IWechsler et al.U2002l) used N-body simulations to study the 
mass accretion histories of dark matter halos observed at var- 
ious redshifts. One of the main results of their work was to 
show that the cumulative mass accretion in individual dark 
matter halos can be well-described by a simple function: 



M(a) = M exp 



-2a, 



(t-0 



(1) 



where a = (1 +z) _1 is the expansion factor and Mq is the mass 
of the halo as observed at the epoch «o = (1 +zo)~ l ■ Given a 
halo mass and redshift, this function has only one parameter, 
a c , which is a characteristic formation epoch for the halo. The 
typical value of a c varies as a function of halo mass, with more 
massive halos having earlier formation times. At a fixed Mo, 
the distribution of a c is roughly log-normal, with a scatter of 
Alog 10 a c = 0.14. Of note is the interesting fact that the a c 
distribution is independent of oq for halos of a fixed mass. 

Wechsler et al. (2002) showed that the median value of 
the formation epoch parameter, a c , can be determined for a 



halo of mass Mq using a simple model based on spherical col- 
lapse. Specifically, a c is set by the epoch at whi ch M*, the 
typical collapsing mass (e.g. lLacev & Colelll993l) . equals a 
fixed fraction F = 0.015 of the halo mass Mq: 



M*(a c ) = FMq. 



(2) 



As expected, this expression for a c is independent of oq. A 
Milky-Way type halo, for example, has a median a c ~ 0.30, 
corresponding to a look-back time of ~ 10.5 Gyr, which is 
consistent with the expected epoch of peak accretion dis- 
cussed above. 

3. STAR FORMATION AND EXAMPLE HISTORIES 

We will use three example mass accretion histories to ex- 
plore expected differences between halo stars and surviving 
satellites. Since each system under consideration is a low- 
mass object with efficient gas cooling, we assume the cold 
gas inflow rate, g(t), tracks the dark matter accretion rate: 



g(i)=f m h(t-t ia ). 



(3) 



Here, / gas is the fraction of the dark matter halo mass in the 
form of cold, accreting gas, and h(t) = dM/dt is the halo 
growth rate determined by taking the derivative of Equation 
(1) with respect to cosmic time. The time lag t-t m between 
h(t ) and g(t ) accounts for the finite time it takes for gas to set- 
tle into the central galaxy after being accreted onto the halo. 
We assume this settling occurs in roughly a dynamical time 
at the virial radius: t m = /?h/v C irc — 2Gyr (1 +z)~ 3//2 . We set 
/gas = 0.02, in acc ord with the cold baryonic mass fraction in 
observed galaxies jBell et alJ2003l) . 

The star formation rate, tp(t), in each system is set using the 
simple relation 



M„ 



(4) 



where M gas is the total cold gas mass in the system and the 
star formation timescale, r*, is a free parameter. In the mod- 
els we explore below we assume for simplicity that star for- 
mation is truncated soon after each satellite halo is accreted 
onto the Milky Way host as a result of ram-pressure strip- 
ping from the backgroun d hot gas halo ( Mo ore & Davisl!*)^ 
Mailer & Bullock 2004). Of course, this assumption is over- 
simplified and cannot be accurate in detail, but it allows us 
to capture the general differences we expect from one exam- 
ple to the next without addressing the complicated interplay 
between orbital evolution and star formation. Ram-pressure 
will be more effective at stripping gas from low-mass objects 
with shallow potential wells (dSph) while realistically dlrr- 
size halos will tend to retain their gas more efficiently. Ad- 
ditionally, for a given mass, typical encounter velocities are 
higher at higher redshifts. Surviving dwarf satellites will un- 
dergo correspondingly less ram-pressure stripping than their 
counterparts that are accreted into the halo early and subse- 
quently disrupted. However, for massive systems that were 
accreted early (a "halo-progenitor", for instance) our assump- 
tion of star formation truncation will not qualitatively affect 
our results compared to a model where star formation is al- 
lowed to continue until the galaxy is tidally destroyed. Early- 
accreted systems are typically destroyed very soon after they 
merge (over a ~ 1 Gyr timescale, a typical dynamical time for 
a z~l halo). Nevertheless, we acknowledge our assumption 
of truncated star formation is simplified. We aim to explore 
the implications of varying models of star formation within 
satellite halos in subsequent work. 
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The gas mass in Equation[4]will not simply be the cumula- 
tive integral of gas inflow rate, g(t ), as we allow for the effect 
of star formation feedback and associated gas loss within each 
satellite halo. The system of equations we use to track star 
formation, gas loss, and metal enrichment is presented in §4. 
Including this gas outflow, we find that the choice r* = 6.75 
Gyr allows us to match the observed mass-to-light ratios in 
dlrr-sized objects, and, when incorporated into our full mass 
accretion history models of the Milky Way, the same choice 
gives a good match to the total Mi lky Way stellar halo mass 
as well as its radial profile (Bullock & Johnston 2005). The 
inclusion of feedback in our model helps regulate the star for- 
mation rate in dwarf galaxies. As a result, we must adopt a 
smaller star formation timesca le than that require d in the sim- 
ilar model utilized by Bullock & Johnston ( 2005) which does 
not explicitly use a prescription for feedback from star forma- 
tion. Before moving on to discuss our treatment of outflow 
and metal enrichment in detail, we first introduce our three 
example mass accretion histories. 

Dwarf Irregular Satellite 

For our surviving, "dlrr" satellite halo example, we use a 
dark matter halo of virial mass of Mo = 6 x 1O IO M , which 
we assume is accreted into the Milky Way halo at ao = 0.79. 
This accretion epoch corresponds to a look-back time of 3.1 
Gyr, a typical accretion time for surviving, massive satel- 
lites (A. Zenter & J. Bullock, private communication) and 
close to that expected for the LMC from gravitational and hy- 
drodynamical simulati ons of its orbital evolutio n within the 
Milky Way potential (Mastro pietro et al J 120041) . A halo of 
this mass at ao = 0.79 has a median maximum circular ve- 
locity of Vcirc = 72kms _1 (Bullock et al. 2001), but we expect 
that dark matter mass loss will leave this system with a circu- 
lar velocity closer to that of the SMC or LMC (~ 60km s" 1 ; 
Ivan der M arel et aT]|2002|) after three billion years of evolu- 
tion (iBullock & Johnstorl2005l) . 

We assume that the dark matter mass accretion history 
of our dlrr halo follows the form of Equation 1, with a 
characteristic formation epoch of a c = 0.231 (Equation 2), 
corresponding to a look-back time of 11.6 Gyr. Note that 
the time between "formation" at a c and its accretion into 
the Milky Way dark halo at ao, this system will have had 
~ 8.5 Gyr to form stars - ample time for both Type la and 
Type II enrichment. After applying the star formation law 
discussed above and the feedback prescription described in 
§4, the final stellar mass and gas mass in this system are 
M* = 3.9 x 1O 8 M and M gas = 4.7 x 1O 8 M , respectively. 
Allowing for a small amount of g as loss from ram pressur e 
stripping after accretion (see e.g. Mastropietro et al. 2004), 
these numbers are in line with gas mass an d stellar mass 
measurements for the SMC and LMC (e.g. iGrebell 1200 It 
iGrebel etal.12001. 

Dwarf Spheroidal Satellite 

Our Dwarf Spheroidal example consists of a relatively low- 
mass halo, Mo = 5.6 x 1O 8 M , accreted at ao = 0.667 or at a 
look-back time of 5 Gyr. The accretion time is typical for that 
of surviving satellite halos of this size (A. Zentner & J. Bul- 
lock, private communication). The halo's formation epoch is 
a c = 0.08, corresponding to a look-back time of 1 3 . 1 Gyr. This 
system has a maximum circular velocity of v c ; rc = 20km s" 1 
at accretion, which will likely be reduced to v C ir C ~ 15kms _1 
afte r five billion years of evolu tion within the galaxy poten- 
tial (B ullock & Johnstonl2 005). This circular velocity agrees 



with circular velocity estimates for dSph s atellite galaxies in 
the Local Group (e.g. IZentner & Bullockl 120031) . The time 
available for star formation between a c and ao is ~ 8 Gyr for 
this system. 

It is now well-known that if ACDM is correct, only roughly 
one-in-ten of the dwarf-size satellite hal os can host a dwarf 
galaxy, with the rest remaining dark (iMoore et alJ 1999; 
Klvpi rTet alJI 1999ft . A natural solution to this problem relies 
on the fact that halos of this size (v c i lc < 30km s" 1 ) are un- 
able to accrete gas after the universe becomes reionized (e.g. 
IBullock et alJ200 0). We account for this effect by setting the 
gas accretion rate g(t) = 0.0 after the epoch of reionization, 
a re . The value of a c we have adopted for our dSph example 
is ~ 2cr earlier than the median expected for a halo of its size 
(via Equation 2). We have made this choice in order to ensure 
that this halo collapses before cosmic reionization and is able 
to host a visible dwarf galaxy. Our assumption about the ef- 
fect of reionization on dwarf galaxy formation allows our pic- 
ture to self-consistently address the "dwarf satellite problem" 
without modifying the tenets of ACDM. For concreteness, we 
choose a le = 0.09, but our results are insensitive to this choice. 

The stellar mass and gas mass within this system at 
accretion is M* = 1.1 x 10 6 M Q and M gas = 3.8 x 10 5 M Q 
respectively. The stellar mass is consistent with that of 
dSph galaxies in the local group, and its gas mass fraction 
is lower than that in the dlrr example, also in acco rdance 
with observed trends (e.g. lMateolll998t IGrebel et all 120031 
and references therein). We expect that the gas mass in this 
system will be further reduced after accretion as a result of 
efficient ram pressure stripping from the shallow host halo 
potential well. 

Stellar Halo Progenitor 

As argued in §2.1, an d discussed in more detail in 
IBullock & Johnston! ( 120051) . the hierarchical scenario predicts 
that most of the mass in the stellar halo was contributed by 
several massive subhalos that were accreted onto the Milky 
Way host long ago. We choose to model one such sub- 
halo as our example "stellar halo progenitor" a dark halo 
with a virial mass equal to that of our dlrr example above, 
M = 6 x 1O 1O M , but which was accreted 9 Gyr in the past, 
at ao = 0.417. Because this subhalo has the same mass as 
the dlrr example, the stellar halo progenitor has the same for- 
mation epoch, a c = 0.231. However, because the subhalo was 
accreted (and destroyed) at a much earlier time, the time avail- 
able for star formation in this system is only ~ 2.6 Gyr, and 
we therefore expect that the ~ 1 Gyr required for Type la en- 
richment will limit significant Type la enrichment in this sys- 
tem. 

When we apply our star formation and feedback model, the 
stellar mass at accretion in this system is M t = 2.4 x 10 8 M Q 
and the cold gas mass is M gas = 7.2 x 1O 8 M . The high gas 
fraction is a result of have a shorter time available to turn gas 
into stars. Indeed, since a majority of the baryons in the Milky 
Way system were accreted in early events of this size (§2.1), 
gas-rich mergers of this type are likely fundamental in fuel- 
ing the growing Milky Way bulge and disk components. Al- 
though a large number of lower-mass subhalo accretion events 
will likely contribute significant substructure to the Galactic 
halo, we expect that ^4-5 mergers of subhalos similar to our 
halo progenitor model are responsible for most of the mass of 
the Milky Way stellar halo. 

In Figure ^ we plot star formation histories as a function 
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FIG. 1 . — Model star formation histories of our halo-progenitor example 
(solid line), dlrr (dotted line), and low-mass dSph galaxy (dashed line). The 
halo-progenitor model forms stars quickly until it is accreted into the stellar 
halo ~ 9 Gyr ago. The surviving dlrr galaxy forms stars more gradually, until 
its recent accretion into the galaxy ~ 3 Gyr ago. The dwarf spheroidal galaxy 
has its gas infall suppressed by reionization early-on and converts its gas mass 
into stars until its accretion into the galaxy ~ 5 Gyr ago. Both the dlrr and 
dSph are expected to survive as Milky Way satellites until the present day. 



of look-back time for each example history: halo-progenitor 
(solid line), dwarf irregular (dotted line), and dwarf spheroidal 
(dashed line). The break in the dSph track at ~ 13 Gyr is as- 
sociated with suppression of gas infall at reionization. From 
Figure \l\ we see that the halo-progenitor star formation his- 
tory is truncated and more rapid compared with the surviving 
dlrr. The history of the dSph example is intermediate between 
the two. We now address how these differences will manifest 
themselves as chemical signatures in the different stellar sys- 
tems. 

4. STELLAR YIELDS 

In order to track the time-dependent chemical content of 
dwarf galaxies during their formation, we model the chemical 
yields from stellar populations, including Type la supernovae, 
Type II supernovae, and stellar winds from intermediate mass 
stars. The rate of enrichment is calculated analytically and 
the yields for individual chemicals are selected from the liter- 
ature. 

4. 1 . Initial Mass Function 
We define our initial mass function (IMF) as 



dN 

4>(m)= — =B{m/M Q r 
am 



(5) 



which provid es the differential num ber of stars of mass m. 
We adopt the iKrouoal d200Tl 12001 IMF, which has a mass- 
dependent slope s given by 



0.01 < m < 0.08 
0.08 <m <0.5 
m > 0.5 

We set the constant B to normalize the IMF to give 




m<p(m)dm = 1M, 



(6) 



(7) 



as we consider stars of mass 0.1 <m< 100. Our normal- 
ization provides <fc(m) in units of Mq and allows us to calcu- 
late our yields per unit solar mass of stars. We note that for 
ste llar masses m > 0.5M© the Kroupa IMF is very similar to 
the Salpeter ( 1955) IMF that has a constant slope of s = 2.35, 
though, as discussed by Kroupa (2002), uncertainties in de- 
termining the IMF of massive stars may mask a steeper slope 
of s w 2.7 for stars above 1 M©. 

4.2. Type II Supernovae 

Given an IMF, we can track the cumulative number of Type 
II supernovae that have occurred per solar mass of star forma- 
tion as a function of time t as 



Mi. 



N n « t) = 



<j)(m)6m. 



(8) 



MM>Mn,^ 

Here, M.*{t) is the mass of a star with stellar lifetime t and 
M-n.mm is the minimum mass star producing a Type II su- 
pernova and yW//,max is the maximum stellar mass that will 
produce Type II yields. We assume that all stars more mas- 
sive than 7W//,max = 50M© will collapse completely into black 
holes and we use A4// m i n = 10M©. 

For the stellar lifetime as a function of star m ass, we adopt a 
dual power- law fit to the Schall er etal] fl992) stellar lifetimes 
of the form 

f * ( - M)=10 - 5 (^) 3 + -° 35 (^) 55 ^ 
with t* in Gyr. The the inverse quantity M.±{t) and its deriva- 
tive dA/U/df* can be obtained from this relation. The rate of 
Type II supernovae per solar mass in each stellar population 

is 

KiK0 = -^flMt(0], (10) 

df* 

with no further normalization being necessary as all massive 
stars between Mu. mm and Ain. max are considered to produce 
supernovae. 

The instantaneous chemical yield of an element X from 
Type II supernovae as a function of time after a stellar pop- 
ulation forms will then follow 

Y XM (t) = y XJ ,[M(t)]R„(t), (11) 

where y xl \M.] is the Type II yield of element X from a star 
of mass M. In prac tice, we use the Type II y ields from 
Tsuiimoto et al. ( 1995) and Thielemann et al. ( 1996), interpo- 
lating between their simulated models in order to obtain yields 
as a continuous function of mass. 

Given the above model, Type II supernovae provide enrich- 
ment from f*(50) = 4.1 Myr until f*(10) = 20.4 Myr after each 
stellar population forms. 

4.3. Type la Supernovae 

In calculating the rate of Type la supernova , we fo llow the 
formalism developed by iGreg gio & Re nzinil Jl983h . which 
postulates that the progenitors of Type la supernovae are bi- 
nary star systems. The rate of Type la supernovae is set by the 
evolution of the secondary star in each binary system. 

We define the number of Type la supernovae that have oc- 
curred before time t as 



0.1 



Mhfim 


' 1/2 




N,«t)=A J 4>{M h ) 


J 2 1 ^(l+ 7 )/Al/i 




M b . mt 


Mini' 
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where A is a normalization parameter. The outer integral sums 
over the total binary mass, Mb = M\ + A^2,with Mi < M\, 
and the term in brackets describes the distribution of sec- 
ondary mass fractions /i = Mi/Mb- The exponent 7 = 2 
sets the spectrum of secondary masses. The upper limit on 
the total binary mass is set by twice the maximum primary 
mass that allows the development of a degenerate C/O core 
before filling its Roche lobe, Mb.max = 16M S . The lower 
limit forces M 2 < Mi via Mh inf = MAX\Mo(t) Mh m m] with 
Mb.min = 3M© following Greggio & Renzinjj (119831) . As in 
§3.2, Mi{t) is the mass corresponding to the stellar life- 
time t. The lower limit on the secondary mass fraction is 
/Lti„f = MAX[7W 2 (0/^, i.M b - jM fejmax )/Mfe]. In o rder to ap- 
proxim ate local supe rnovae rates as determi ned by Tammann 
(119821) . we follow the lMatteucci & Greggiol J 19861) Type I SN 
normalization. Adjusting the Matteucci & Greggio ( 1986) re- 
sult for our chosen IMF, we find A w 0.0956. 

Given the cumulative number of Type la supernovae as a 
function of time f , the corresponding rate is simply 



R,(t) = 



dNi(< t) 
df ' 



(13) 



The instantaneous chemical yield for each element X from 
Type la supernovae is then Y XI (t) = y x ,Ri(t). For the individual 
yield per supernovae yxj, we adopt the yield s of the Nom oto 
W7 Type la supernova model (Nomoto et al. 19841 11997ft . In 
this model the Type la supernovae provide enrichment for t > 
f+(8) = 31.7 Myr after each stellar population forms, with the 
peak of the Type I rate occurring at t ~ 1 Gyr. 

4.4. Stellar Winds 

We follow the chemical enrichment from stellar winds by 
tracking the death rate of stars and assuming the majority of 
the stellar mass loss occurs during asymptotic giant branch 
thermal pulsing near the end of the stellar lifetime. Our effec- 
tive rate of stellar winds is then 

AM 

Rw(t) = --7-<KM), (14) 
df* 

which, like the Type II supernovae rate, is set by the IMF 
and stellar lifetimes. The instantaneous chemical yield of 
an element X f rom stellar winds is then Y vu ,(t ) = y VI , ,Rw(t)- 
We utilize the van den Hoek & Groenewegen (1997) time- 
dependent stellar wind yields yx,w, interpolating over mass 
and metallicity for stellar masses between 0.8 -8Mq. Stel- 
lar winds provide enrichment in this model continuously after 
f*(8) = 31.7 Myr. 

5. CHEMICAL ENRICHMENT AND FEEDBACK 

5.1. Tracking Abundances 

We aim to track the abundances of H, He, Fe, and the a- 
elements O and Mg using the three example accretion histo- 
ries discussed in §2. To do so, we follow the evolution of 
stellar mass M*, gas mass M gas , and mass M x in each element 
X using the system of equations: 

dM* 

*=m-Y(t) (15) 



dt 

dM™, 

—J^=-^(t) + Y(t)+g(t)-e(t) 
at 

dM v 

-f x ^(t) + Y x (t) + h x g(t)-f x e x (t). 



dt 



(16) 



(17) 



rate, g(t), tracks the halo mass accretion rate. The quantity 
Y(t) is the total rate of mass input into the gas from super- 
novae of both types and stellar winds, summed over all stellar 
populations. Similarly, Y x (t) is the instantaneous input rate 
of the specific element X, again summed over supernovae and 
winds for all stellar populations. The gas mass fraction in each 
element X is f x = M x /M gas , and h x is the elemental mass frac- 
tion in newly accreted gas. The ejection rate for each element 
is chosen to follow the star formation rate as e x (t) = w x ijj(t), 
with w x representing the efficiency of gas blow-out for ele- 
ment X. The total gas ejection rate in the second equation is 
then the sum over all elements X, e(t) = Y.f x e x (t). 

Each star formation history begins by accreting low- metal- 
licity, alpha-enriched gas, with [Fe/H]=-4.0 and [a/Fe]=0.3. 
Throughout the history of each galaxy, we assume that the 
infalling gas retains this initial Type II-enriched abundance 
pattern (this sets h x in Equation 17), although, in detail, our 
model results are insensitive to the abundance pattern of the 
low-metallicity infalling gas. In the results presented in |6j 
we define our abundanc es relative to t he solar abundances 
measured by Grevesse & Sauval ( 1998). We turn now to a 
more in-depth discussion of our treatment of blow-out feed- 
back. Our feedback prescription is important for explaining 
the abundance pattern in low-mass, dSph-type galaxies. 

5.2. Feedback and the Mass-Metallicity 
and Mass-v c i vc Relations 

Gas and metal ejection from supernovae winds almost 
certainly plays an important role in shaping the properties 
of ga l axies - especially dw arf-sized galaxies llDekel & SiUd 
1986). lDekel & Woolll2003l) (DW) have used a compilation of 
dwarf galaxy data to show that the stellar masses and metallic - 

2/5 

ities of dwarf galaxies follow a fairly tight relation, Z oc MJ , 
and have argued that blow-out feedback can naturally explain 
this trend. Although si milar r elations were n oted previously 
(e.g.. iBrodie & Huchralfl99l IZaritskv et alJll994l) . we refer 
to this trend as the "Dekel-Woo relation", and use it to mo- 
tivate and normalize our prescription for galactic supernovae 
ejecta. The feedback prescription we present below in fact re- 
produces the Dekel-Woo relation, but also gives rise to and ex- 
plains the [a/Fe] abundance pattern trends observed for low- 
mass dwarf galaxies. We note that our prescription for gas and 
metal ejection does not strongly affect our conclusions regard- 
ing the stellar halo progenitor abundance pattern nor those of 
our surviving dlrr example. The robustness of the model pre- 
dictions for these subhalos results because the dark halos of 
these systems are relatively massive. However, ejection plays 
an important role in the star formation and enrichment history 
of our lower-mass dSph example, where the typical energy of 
supernovae ejecta is significant compared to the binding en- 
ergy of the halo. 

We can estimate the efficiency of gas and metallicity ejec- 
tion from dwarf halos, w, by taking the ratio of specific energy 
of gas that escapes the central dwarf galaxy as a wind, £ w , to 
the specific binding energy of the halo: 



f 



2 cue 



2 circ 



(18) 



As discussed in §2.3, the star formation rate is assumed to 
track the mass in cold gas as tp(t) = M^/t^, and the gas infall 



where v c i, c is the circular velocity of the halo in question. The 
energy of the gas that leaves the central galaxy will depend 
on A^sn^ the number of supernovae per unit mass in the stellar 
population and Z?sn> the energy per supernovae. The factor 
r] is the fraction of supernovae energy converted into kinetic 
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energy, and e is the fraction of the input kinetic energy that 
is retained by the gas that actually escapes the galaxy as a 
wind. The kinetic energy fraction 77 will depend on the speed 
at which supernovae bubbles overlap compared to the radia- 
tive timescale o f the gas in the gal axy, and therefore on the 
supernovae rate (Dekel & Silk 1986). It typically will not ex- 
ceed a few percent, and is commonly taken to b e 5% in cos- 
mological simulations of galaxy formation (c.f. lAbadi et'al] 
2003). With this, typical numbers give w w 22er] 5 v^Q, where 
■q 5 = ?7/0.05 and V50 = v c i rc /50 km s" 1 , and we have used 
A^sn ~ 1.1 x 10" 2 Mq, which follows from our chosen IMF 
and stellar population modelling. 

The fraction of kinetic energy e retained by gas that escapes 
as a wind will vary between elements. While H and He will 
primarily not be entrained in the wind because of interactions 
with the surrounding ambient medium, e <C 1, metals will be 
preferentially blown -out of galaxies in the metal-rich super- 
novae ejecta, e ~ 1 JMac Low & Ferrara 19991) . Furthermore, 
the rate of Type la supernovae per stellar population is quite 
slow, dragging out over Gyr timescales, while the Type II su- 
pernovae rate is rapid. We therefore expect that the bubble 
overlap time and kinetic energy fraction to be slightly lower 
for Type la compared to Type II ejecta: r] a < r\p e . Motivated 
by this physical picture, we set the gas blow-out efficiency of 
gas, Fe, and a-elements in halos of varying circular velocities 
to be 

- g as(v circ ) = 0.15x(^ zr )" (19) 

^ ( ^ )=14 - 4x (50k^)" (20) 

-(vc ir c)=16.8x(^ rr )" 2 . (21) 

We have adopted the standard kinetic energy fraction 77 = 0.05 
for gas and Fe, but a slightly lower value, 77 = 0.035, for the a 
elements. We have used e = 0.01 for the gas and e = 0.8 for 
the metals. 

While our choices for Wf e and vv ga s are of roughly the ex- 
pected size, we have chosen their precise values in order to 
reproduce the observed - v c \ rc relation for dwarf gal ax- 
ies in the local group (DW; see Bullock & Johnston 2005) as 
well as to approximate the Dekel- Woo mass-metallicity re- 
lation. Figure [2] shows the stellar-mass metallicity relation 
estimated by DW as the solid line, with the dashed lines in- 
dicating the I— a scatter in metallicities about the relation. 
The two solid circles and the open circle show our model re- 
sults for the dSph, dlrr, and stellar halo progenitor examples, 
respectively. The error bars reflect the 68% spread in metallic- 
ity range for stars in each model galaxy. In order to more fully 
compare with the observed relation, we have also included re- 
su lts for a range of surviving satellite halos from the models 
of Bullock & Johnston (2005). The mass accretion histories 
of these satellites were generated in the same way discussed 
§2, except now we span a more complete range in v c ; rc . The 
model galaxies generally reproduce the observed trend, with 
good agreement considering the uncertainty in the data and 
the difficulty in assigning precise stellar masses to these ob- 
served systems. We retain the stellar feedback calibration that 
provides this agreement for the proceeding calculations. 

Note that there is very little freedom in our parameter 
choices. Recall that f* and / gas in our mass accretion model 
are already fixed to reproduce the stellar mass of the Galactic 
halo and the cold baryon fraction in galaxies. Once these are 
given, w is constrained by the stellar mass-v c ; lc relation, and 



0.0 
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FIG. 2. — The stellar mass - me tallicity re lation for local dwarf galaxies. 
We plot the observed trend from Dekel & Woo 1 2003) (solid line), the region 
containing 70% of the observed data (dashed lines), and the calculated rela- 
tion for a range in mass of model dwarf galaxies (crosses). The error bars 
on the model calculations delineate the ~ 1 — a metallicity range (68%) of 
stars in each system. We use this calibration to constrain the strength of stel- 
lar feedback in our model. We also plot our models for surviving dSph and 
dlrr galaxies (filled circles) and the halo-progenitor model (open circle) for 
comparison. 



wpe by the approximate stellar mass-metallicity relation. Fi- 
nally, the parameter w a is chosen to reproduce observations 
of the chemical abundance pattern of dSph type galaxies once 
the above list of complimentary observational data is invoked 
to constrain our other free parameters. In the next section 
we show that our choices for the wind efficiencies indeed re- 
produce observations of the chemical abundance pattern of 
dSph-type galaxies. We argue below that the chemical con- 
tent of low-mass dwarf galaxies should then show trends as 
a function of stellar mass and circular velocity. Specifically, 
we expect that high-mass systems will tend to have near solar 
abundance patterns, while low-mass dwarfs will have inter- 
mediate [a/Fe] ratios. 

6. RESULTS 

The three panels in Figure [3] show the [a/Fe] vs. [Fe/H] 
tracks computed for our stellar halo progenitor (top), dSph 
example (middle), and dlrr example (bottom). The tracks cor- 
respond to the star formation histories presented in Figure ^ 
In all cases, solid lines show [O/Fe] and dotted lines show 
[Mg/Fe]. The final stellar masses (at accretion) for the sys- 
tems are indicated in the bottom left of each panel. The thick 
region along each line corresponds to the abundance region 
inhabited by the middle 68% of stars in each system (the 68% 
spread by mass about the median). 

In the upper panel, our massive (Mo ~5x 10 I0 M©), early- 
forming, stellar halo progenitor forms stars quickly until it is 
accreted into the Galactic halo ~ 9 Gyr ago. Its gas enriches 
to [Fe/H]= -0.7 and an a-abundance of [a/Fe]=0.1 just be- 
fore accretion, while its stars typically have a metallicity of 
[Fe/H]~ -0.9 and [a/Fe] ~ 0.2. Upon accretion, we expect 
this system to quickly disrupt, and contribute its stars to the 
stellar halo (§2). Over the predicted chemical evolution tracks 
we plot abunda nce data (cros ses) for galactic halo stars from 
Table 2 of IVenn et al.l J2004I) . from which we select all stars 
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-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 

[Fe/H] 

FIG. 3 . — The [a/Fe] vs. [Fe/H] abundance of our halo-progenitor (top panel), dwarf spheroidal (middle panel), and dlrr-type galaxy (lower panel). Top: The 
halo-progenitor model retains a super-solar [a/Fe] abundance pattern, resembling the stellar halo data (crosses, from Venn et al. 2004). Middle: The abundance 
pattern of the dwarf spheroidal galaxy is set by the strength of the stellar feedback, which forces its [a/Fe] abundance below that of halo stars. For comparison, 
we overplot abundance data of similarly-sized dSph galaxies from Venn et al. 1 2004) as crosses. Bottom: The majority of the star formation in the massive dlrr 
occurs at near-solar values of [a/Fe] and at higher metalhcities, similar to stellar abundances in the SMC (stars and crosses from lVenrl 19991 IHill et all 1993. 
respectively). Note that we expect extremely metal poor stars in the SMC and LMC to have enhanced, halo-like [a/Fe] ratios. 



with Galactic halo membership probabilities greater than 50% 
5 . The halo-progenitor dwarf galaxy abundance pattern is con- 
sistent with the overall abundance pattern of the Galactic halo. 
Agreement between our model and the data arises because 
the stellar halo progenitor is accreted and destroyed before 
significant Type I enrichment can occur. More generally, we 
expect that the super-solar a-abundance of halo stars is inher- 
ited naturally from the star formation and chemical enrich- 
me nt histories of these massive, ea rly-accreted dwarf galax- 
ies iFont et alJ l2005. in preparation). As discussed in §2, the 
- 1.5 x 10 9 M Q stellar mass of the halo JCarnev et alii 19901) 

The data was ori ginally obtai n ed bv IFulbrighH 120001 12001): 
IGratton & Snedenl ROM (l991t H994l: IHanson et all 11998ft: llvans et alJ 
120031): McWilliam eTall l\993l : IMcWilliaml I199ST : iRvan et alJ 1199^K 

Stephens & Boesg aann2002l) . 



is expected to be dom inated b y stars fr om disrupted systems 
of this type jBullock & Johnstonl20 05). 

In the middle panel of Figure |5J our low-mass (Mo ~ 
5 x 10 8 M Q ), dSph example forms stars more gradually, and 
is accreted onto the Milky Way halo ~ 5 Gyr ago. Metal en- 
richment is suppressed in this system compared to the other 
two panels because the dSph halo is more loosely bound and 
supernovae winds can more easily escape (§4.2). Since the 
overall ejection rate is higher, the relative efficiency of blow- 
out for Type I and Type II supernovae ejecta is more appar- 
ent and creates the break in the trajectory in the dSph track. 
For comparison , the cros ses show a-abundance data com- 
piled by Venn et al. (2004) for Local Group dwarf spheroidal 
galaxies. The data reflects the well-established disparity be- 
tween the [a/Fe]~0.3 abundance ratios of Galactic stellar 
halo stars and the [a/Fe]^0.0 - 0.2 abundance ratios of lo- 
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cal, low-mass dwarf galaxies at similar metallicities. It is en- 
couraging that our model dSph evolves on an a-abundance 
track that overlaps with data of observed dwarf spheroidal 
galaxies, with its final stellar abundance peaking at [a/Fe]~ 
0.05 and [Fe/H]~ -1.65. We emphasize again that this agree- 
ment stems from supernovae blow-out. The differences be- 
tween the a-abundance pattern for our halo progenitor and 
our dSph example primarily results from the relative poten- 
tial well depths of their host dark matter halos. However, we 
note that differences between the a-abundance pattern of spe- 
cific dSph galaxies may reflect different efficiencies for blow- 
out of Ty pe la and Type II supernovae ejecta between sys- 
tems (c.f. Lanfranchi & Matteucci 2004) than what we adopt 
above. In principle, differences in blow-out efficiencies be- 
tween systems will depend on the complicated hydrodynam- 
ical problem of following the overlap of supernovae bubbles 
and cannot be properly addressed by our simple modelling. 
Similarly, our choice of initial mass function affects the num- 
ber of supernovae per unit mass of stars formed, but such al- 
terations to our modeling can be accomodated through our 
feedback calibrati on and similar results can be obtained for 
e.g. the Salpeter (1955) IMF. We therefore rely on the ob- 
served trends of a-abundance in dSph galaxies, the Dekel- 
Woo relation, the total stellar mass of the Galactic halo, and 
the stellar mass-v c j lc relation expected for Local Group dwarfs 
to constrain our singular choices for the blow-out efficiencies 
that we utilize in our modelling. 

Finally, the [a/Fe] vs [Fe/H] track for our massive (Mo ~ 
5 x 1O 1O M0) dlrr example evolves in a similar manner to that 
of our stellar halo progenitor case, except that it is accreted 
much later and can continue evolving for several more billion 
years. Because of this longer, slower star formation history, 
the effect of Type I supernovae is more pronounced than in the 
halo-progenitor case. The gas metallicity enriches to [Fe/H] 
= -0.6 with an a-abundance of roughly solar, while the in- 
ner 68% of stars by mass span a metallicity of [Fe/H] = -0.95 
to [Fe/H] = -0.6 and range from [a/Fe] = 0.2 to near solar 
values. Indeed, it has abundance patterns typical of massive 
satellite galaxies in the loca l group, as illustrated by the over- 
pl otted from LM C stars bv lHill et all i 19951) and SMC stars 
bv iVennl {1999). Note that we predict that the lowest metal- 
licity stars in galaxies like the SMC and LMC should have 
abundance patterns similar to those in halo stars. The similar- 
ity of a-abundance patterns in dwarf galaxies at early times 
suggests we cannot rule out a scenario in which the Milky 
Way stellar halo was completely comprised from the destruc- 
tion of dSph galaxies before they were significantly enriched 
by Type I supernovae. Ho wever, the statistics o f dark matter 
halo accretion histories (see Zentne r et all2005l) suggest such 
a scenario is unlikely in the context of hierarchical formation 
of the stellar halo. Furthermore, dynamical modelling of the 
formation of the stellar halo might rule out such a scenario 
through differences in the resultant spectrum of phase-space 
correlations. 

Overall, the differences in the shape of the abundance pat- 
tern tracks between the halo, dSph, and dlrr models are 
straightforwardly attributable to the mass-dependent strength 
of the feedback implementation. The stronger winds in the 
low-mass dSph model cause the [a/Fe] abundance to drop 
near [Fe/H] = -2 as ^70% of the metals produced by star for- 
mation have been expelled from the system by this time. At 
the same metallicity, the halo progenitor and dlrr models have 
only lost ^7% and ^9% of their metals, respectively. For the 
dSph model, Type I SN begin to dominate over Type II SN by 



this metallicity and the retention of Fe relative to a-elements 
correspondingly increases. The relative strength of the a- and 
Fe-wind effic iencies will vary between indiv idual dSph galax- 
ies (see, e.g.. lLanfranchi & Matteucci 2004) to produce their 
various abundance patterns, but the [a/Fe] tracks of the more 
massive halo progenitor and dlrr systems will be fairly in- 
sensitive to these details as their less efficient winds begin to 
strongly change their abundance patterns at characteristically 
higher metallicities. 

7. CONCLUSIONS 

By combining cosmologically-motivated star formation and 
gas accretion histories for dwarf galaxies with modelling of 
the chemical yields from stellar populations, we demonstrate 
that the abundance patterns of the Galactic halo, low-mass 
dwarf spheroidals, and more massive dwarf irregular galax- 
ies can be understood within the context of cosmologically- 
motivated models of the stellar halo and dwarf galaxy forma- 
tion. A straightforward prediction of ACDM is that most of 
the stars in the stellar halo were formed in dlrr-sized dark mat- 
ter halos that were accreted at early times and subsequently 
destroyed. The destroyed, proto-halo systems are short-lived, 
form their stars rapidly, and are enriched predominately by 
Type II supernovae. They are disrupted before the iron input 
from Type I supernovae can drive them to higher metallicities, 
leaving a chemical abundance pattern that matches observa- 
tions of low-metallicity, a-enriched halo stars. 

Longer-lived, dlrr-type galaxies that survive to the present 
day, form their stars more gradually. These systems undergo 
both Type I and Type II supernovae enrichment. We presented 
a "typical" dlrr-type accretion history, which enriched to a 
moderate metallicity with an abundance pattern similar to that 
of the SMC and LMC. We predict that lower-metallicity stars 
in the SMC or LMC will likely show similar abundance pat- 
terns to those of halo stars. 

Low-mass dSph-type galaxies are likely systems whose 
chemical abundance patterns have been altered significantly 
by supernovae winds, as their shallow potential wells cannot 
retain outflows. Adopting simple assumptions about the ejec- 
tion wind efficiency, we showed that the expected population 
of surviving dwarf satellite galaxies in ACDM galactic halos 
match the stellar-mass/metallicity relation observed for local- 
group dwarfs (Dekel & Woo 2003). We then argued that the 
same stellar feedback is essential in explaining the a-element 
abundance patterns in dSph galaxies, and demonstrated this 
idea by presenting an explicit, cosmologically-motived dSph 
model that follows the abundance pattern of the observed 
dwarf spheroidals successfully. Our expectation is that the 
abundance data for dwarf spheroidal galaxies as a function of 
their stellar mass can be matched by similar, cosmologically- 
motivated star formation histories (Font et al., in preparation). 
We predict that [a/Fe] ratios should gradually approach the 
solar value in dwarf galaxies of increasing stellar mass (or ve- 
locity dispersion). 

As discussed, most of the stars in the Galactic stellar 
halo are likely remnants of early-forming dlrr-type galaxies 
that chemically bear little resemblance to present-day dwarf 
spheroidals. At the same time we expect that a significant 
number of dSph-size galaxies have been accreted and de- 
stroyed over the Milky Way's history. Because of their shal- 
low potential well depths, these low-mass dwarf remnants 
should be chemically distinct from most of the background 
population of the stellar halo. Fossil evidence of these accre- 
tion events may be more easily detected if chemical signatures 
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are used to help identify substructure (e.g. lNissen & Schuster! 
[T997tlFontetal.l2005. in pr eparation). We also note that fur- 
ther constraints on our modeling are supplied by observations 
of the age and metallicity distribution of halo stars. While 
the modeling presented here demonstrates that the chemical 
abundance pattern of typical halo progenitors match the ob- 
served [a/Fe] ratio in the Galactic halo, our model must be 
extended account for the full ensemble of halo progenitors 
expected from the time-dependent mass spectrum of accreted 
substructure predicted by the ACDM scenario for structure 
formation. Such an im provement to our modeling will be pre- 
sented in future work (Font et al. 2005, in preparation). 

While we adopt a feedback formulation commonly used to 
model chemical evolution in galaxies, we note other feed- 
back models may be more physically motivated. For in- 
stance, hydrodynamical simulations of galaxy formation may 
include multiph ase descr iptions of the interstellar medium 
jMcKee & OstrikedfT977h that may more naturally produce 
self-consistent star formation histories for dwarf galaxies 
as the microphysics of th e ISM are better captured (e.g. 
Sprin gel & Hern auist 2003). Hydrodynamic simulations can 
also include the effects of supernovae-driven galactic winds 
formulated in a manner similar to what we employ, but the 
detailed calculations of the expansion and cooling of shocks 
allow a reduction in the net number of free param eters (for in- 
stance the need for e is eliminated, see Springel & Hernauist 
120031) . As these forms of feedback affect the star formation 
history and gas content of a given object, they will also af- 
fect galactic chemical evolution. We also note it is possible 
that some dwarf galaxies with spheroidal component s con- 
tain in termediat e-mass su permassive black holes (Bar th et all 
120041: iGreene & Holl2004l) . and hydrodynamical simulations 
of mergers of galaxies containi ng supermass ive black holes 
with thermal feedback driven by Bondi ( 1952)-style accretion 
( Sprin gel et alJl2004tlDi Matteo et all2005l) demonstrate that 
supermassive black holes may have considerable impact on 
the star formation history and gas content of their host sys- 
tems. Feedback from supermassive black holes may then sig- 
nificantly alter the chemical evolution history of galaxies in 
a manner not considered here. Due to the complex hydrody- 



namical and radiative cooling issues posed by supernovae and 
black hole feedback, we suggest self-consistent hydrodynam- 
ical simulations that include the relevant physics described 
above may be necessary to model more completely the ef- 
fects of feedback on the chemical evolution of dwarf galax- 
ies and the Galactic stellar halo. However, to demonstrate 
the feasibility of a cosmologically-motivated explanation for 
the abundance of the stellar halo we believe our compara- 
tively simple model is sufficient. We note that our chemical 
enrich ment modeling could be inco rporated into the calcula- 
tions of Bullock & Johnston (2005) to form a more complete 
framework to model the detailed properties of the Milky Way 
stellar hal o. We plan to p resent this model in a future paper 
(Font et al. 2005, in preparation). 

We close by mentioning that our model does not capture 
the full complexity of the availabl e observations, no tably of 
heavy elements like Ba, Y, and Eu ( iVenn et all2004l) . and fu- 
ture inconsistencies between the scenario for the formation of 
the halo presented here and observations may arise. Future 
modelling could conceivably help draw stronger connections 
between the very detailed stellar observation and cosmolog- 
ical scenarios for the origin of dwarf galaxies and the stellar 
halo. For instance, as full chemical modeling becomes incor- 
porated into cosmological simulations of galaxy formation, 
questions about the chemical abundance of thick disk stars 
and accreted dwarf systems (c.f. lWvs e]|2004j) will find more 
satisfying answers (see e.g.. [Brook et alJl2004 bl We there- 
fore encourage future attempts to coordinate detailed chemi- 
cal modelling with cosmological calculations to increase the 
interface between cosmology and galactic astronomy. 
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